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SIZING lUBE-FIN SPACE RADIATORS 


INTRODUCTION 


The radiator area required to reject a given amount of energy can be 
calculated by direct application of the Stofan-Bolt/.niann radiation law; 


Q=<r<AFT^ . 

Even though this law Is mathematically simple. Its application to radiation sl/ing 
can become complex. To avoid complexity, equation (l) la sometimes applied 
by assessing the effective temperature, T^^^ of the radiator: 


<rcFT^„ 


The effective temperature of the radiator Is assessed on the basis of experience 
and empirical data. This approach Is normally apidled as a result of quick 
needs by project iiorsonnel. However, this approach does have a ”bullt-ln 
capacity to prcxluce large errors. This results from the fourth power relation- 
ship. A small error in the effective temperature Is multiplied several times 
In the resulting area. 

The sensitivity of this error can be determined (pialltatively from the 
previously mentioned Stefan-Boltzmann relationship. The c. tango in the required 
area with an accompanying change in lemitcrature is 


- '*^j ^dT 
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Normalizing these results by substituting equation (l), 


dA . dT 




(0 


R^tl^ th« dA/ A I. the p.rc«.l d-ng. 1. .«• tb.t 

in Umparttar., dT/T. A nnit percent Increw *“ 

In n fcut unit decrenee of the re«ulred eren. An ^ 

effective temperature by 4 percent will undersize the retired area by 16 par- 
ent. This ira significant error, even to be tolerated In preliminary 
However, It should be noted that the percent change In temperatare 
^ *,;Lynnmlc »mper.n.r.. Alrc. «p-Uon (4) i. . 
iSTtlondilp whoro th. dlfleronco botwoon t*o temperature. 
oractlce, where finite differences are encountered, the actual multlpllcat on 
Tr^r ll gr^Ur than four. Consider the example where the actual eff^tlve 
temnarature is 53 2*F (513.2*R), see Appendix A, and an assumed effective 

toml^r.tur.of68U-FI(404»7)/21. 528.5-R. Thl. 1. . percent temperature 

error of 2.98 percent: 


513.2 


2.98% 


The percent error In area Is 


A -A^ 
e o 




1 7 = 12 . 21 % 


.here A_^ I. c are. reaulttn* from T„. and A^ la the area roaultln* from. T,. 

One of the purpoeea heroin ie to preaent a technltpie lor ••alaln^’ radlatera 

which can be defended with riaoroue onatnoorina analyala. ’'J® * 
preaentad In detail. Includlna a .■omputer program to accompllah thebaalc olalng 

Usk which Is compatible with a phase A/ B study effort. 
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SIZING VERSUS DESIGN 


lUdlator slzlnK can be < hara. teri/od by fin ctflcleney (diacnaarf In next 
aectlon). Radiator deaian ia .Mmctorlred by the conniiuratlon to achieve a 
dealred fln .Blolency. The proccdute la to. firat, alre the 
a dealred fin efflelency; aeeond, thla tin cffuloniy la euaranteod by the >.< m - 
n^ratlon from tthleh welitht ran be ealcnlated. B thla weiRht “ 
then a tradeoff haa to be made between alre and weight (or the thernta 

be reduced)# 

The Importont hlet to ret ognl/.e ia that a relationahlp doea exist b«» eon 
radiator area and weluht, depending u|)on the fin effic iency ( FijSa ). 
relative scale has been selected bctvcon one and tuo since anahtlcal 
show approximately a 2 to 1 inverse relationship between area and welj;ht Fot 
examplT. If it was desirable to reduce the area of a tapered tube-fin conflRura- 
tlon by 100 percent, then the weight must Increase by 100 |>erccnt. The extia 
welKbt la manifested in the extra fin mass rcciulred to achieve a greater fin 

efficiency# 
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Hadiator confium ration and weight/ area behavior. 


IJ liJil! 






The desiKn spectrum represented in Figure 1 is bounded by two config- 
urations. The heaviest is a rectangular duct, having (by definition) a fin 

efficiency of 100 percent. However, this rectangular duct will have a relatively 
small area for a given heat load. 

The other end of the spectrum is a tapered tube-Hn. It has been deter- 
mined by rigorous analysis that this configuration has the optimum area to 
weight ratio. Neither of these two configurations, as depicted, is a practical 
consideration; the first is too massive and the second is structurally weak. 
Practical radiators are somewhere between these two extremes, with fln * 
efflclencies between 75 and 90 percent. Generally speaking, from the tapered 
tube-fln configuration, this represents a 50 percent Increase in weight with a 
12 i^rcent decrease in area |1J. Thus, sizing a radiator not only depends upon 
the heat load, but also upon the allowable weight which results from design 
considerations. This report >s primarily concerned with sizing rather than 
design. Design usually occurs in phase C or D. in a primarily design effort, 
and sufficient data exist to select practical fin efficiencies. This allows the 
sizing process to proceed in support of programmatic decision. Thus, there 
are no long delays in specific radiator sizing and weight assessments. 


FIN EFFICIENCIES AND FIN EFFECTIVENESS 


A tube-fin configuration is practical because of strength and rejection 
capability, as previously mentioned. As a result, fin {.erformance is important. 
At least two criteria exist in the literature for assessing fin pierformance: fln 
efficiency and fin effectiveness. Fin effii iencv is defined as the ratio of the 
actual heat rejected to what would be rejected if the entire surface was at the 

r^t fin temperature, T . The environmental effects are accounted for by the 
sink temperature, T : 
s 
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Fin effectiveness, 8 , is defined as the ratio of the actual heat rejected to what 

would be rejected if the entire surface was at the root fin 

mental factors are accounted for by the net heat flux absorbed by the surface 

from the environment: 




Thus, the relationship between the two efficiencies is 111 


■-Ct 


The characteristic behavior of these efficiencies is usually c-haracter- 
ized by a fin profile number which results from rijtid fin analysis. The profile 
numbers are a dimensionless set of characteristic s which are mdicative of 
geometry, material, and local thermal conditions. The characteristic s of 
'm<l.ncy .re niu..r,ted in KlRure 2. Normally, practical r«llalor. have pro- 
file numbers less than 1.0. Thus, fin efficiency Is always greater than 60 
percent. Typical characteristic values of radiator configurations previously 
discussed are illustrated in Figure 2. 

It is lnnK.rtant to rec c>gni/.e that the fin effic iency, as defined, is for a 
Single root-fin tem.ierature. In an actual radiator, the root fin temperature 
decreases in the direc tion of Row. As illustraleci in Figure 2. the loc-al value 
of efficiency is lowest at the inlet cenditions and inc rease m the direction of the 
outlet. Thus, fin efficleiuy cannot be applicnl directly but must be 
over the radiator aria. Surprisingly, a thermal mcxiei employing fin efficieno 

is not readily available. 

The prcHeduro reviewed herein utill/es fin effectiveness as defined in 
Reference 2. The rationale for selec ting this method is its re^y 
Uv to the prellminarv design function. Assumptions are employed which simp 
the problem for easy equation solving ccmj.utcr techniques. View factors and 


U- J 







Figure 2. Fin offioiency and its variation 







Rad-K* a are Inputs which can be evaluated by other programs such as 1X)IIARP 
If accuracy is not extremely important, view factors and Itad-K* a can be 
evaluated by charts in combination with experience. 


The basis for the method einpluyed is the equivalent width relationship, 
employed by Mackey (2| between a rectangular duct and a tube-fin (*onfiguration. 
The procedure is to calculate the required area of a rectangular duct for a given 
heat rejection load. This area is then modified by the equivalent length rela- 
tionship to establish the area reqiiii ed for the equivalent tube fin configuration. 
To demonstrate this tct'hnique, consider the rectangular duct in Figure 3 with 


rejection area A 


The equivalent tube-fin configuration has an area A 


proporti^mal to ( 1. r- 2h ). The ivlationship between these two areas for the 
d n 


same heat rejection tairability is 
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Figure I. K<|ui\:ilenl ivclangular duct for tulM‘-fin. 
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where L la the equivalent recUngular duct width aa deAned by Mackey [2] : 
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(9) 


Mackey notes that the greatest value of is 1.07 which is for the optimum area 

to weight configuration* The smallest value is 1*00* Mackey provides a chart 

for F as a ftinction of profile number. Thus, if the equivalent area of a rec- 
r 

tangular duct is known* then the adjustment can be made to find the area of a 
tube-fin configuration. Before the assessment is made* an evaluation of the 
configuration and environmental conditions must be made* 

This proc-edure is considered to be valid for rectangular and tapered fins* 
V a more complex configuration is involved* the procedure is still valid; how- 
ever* within the computer program* provisions are made to account for the 
temperature gradient between the heat transfer fluid and the radiating surface* 

The computer program for sizing purposes given in Appendix A calculates 

the area required by a rectangular duct* It is then necessary to manipulate this 

value by the equivalent length concept to establish the required area for a tube- 

fin configuration. Combining equations (8) and (9) with Q.__ equal to zero: 

ABS 



A typical fin effot tivuness for low temperature radiators is 70 percent* Usually 

the ratio of 1. / 1. will be approxiinatciv 2* If* F is 1.04* 
n d r 
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Thus, the menlpulatlon required la very simple to arrive at the desired tube- 
fln configuration. 


SPECIAL RADIATOR RELATIONSHIPS 


There is a special case of radiator design of particular Interest that 
arises when the absorbed flux can be assumed to be zero. Under this assump- 
tion the radiator equation simplifies, and several expressions result which can 
serve as a guide In developing a philosophy for particular radiator problems. 

The first of those is the relationship between heat rejection area and 
thermal load: 



The rejection area is that of a rectangular duc t. The equivalent length modifica- 
tion can be applied If a tube-fln c-onflgu ration Is doslred. The Importance of 
equation (12) Is the sensitivity of the radiator inlet and exit temperatures. It 
is not apparent, but for a given Inlet tcmt>ci aturc, the required rejection area 
decreases as the exit temperature Increases. However, as the exit tom|^rature 
increases, the required flow rate through the system also increases. Those 
facts are illustrated In Figure -I. The ordinate scale has been normalized. 


Mass flow rate is normalized to i:i G87 Ib/hr which occurs at an exit 
temperature of 80*F. Area is normalized to 5:»7 ft* (O = 0,70) which occurs 
at 0* F exit temperature. Those data were actually obUlned from the computer 
program of Aivendix A. The effec ts of the heat transfer film coefficient, as a 
result of flow rate, is ac-counted. However, this is an Insensitive consideration 
fbr Reynolds numbers above :J00. The radiator area is sized by the radiator 
thermal resistance. In practical applications, the pump size orpump power may 

not be allowable, within an allcmable mass flow range, however, there Is some 
flexibility in reducing radiator area. 
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RAOIATOII fXIT TEMKRATlifll 


Figure 4, Influence of exit radiator temperature on area 
and flow rate for a tubo-fln configuration. 


The efftvtive tninperature of the I'adiator ia of interest to the engineer 
even though it has little practical value. The effective temperature is defined 
by equation (2). The primary purpose for presenting a rigorous expression 
for this temperature is to demonstrate, to some level, how errors can occur 
by assessing it by experience or average values: 



The effective temperature as computed by the radiator program is presented in 
Figure S. 
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0 s so 7f 100 

RADIATOR EXIT TEMPERATURE (° PI 


Figure 5, Sensitivity of effective temperature to exit 
temperature and absorbed flux* 

Equation (13) is specially for zero absorbed flux. However, Figure 5 
hao additional data to illustrate how view factor and absorber flux can attect the 
effective temperature. To illustrate ftirthcr the sonsittvtty of effective tem- 
perature, the dashed lino is for an absorbed flux of 30 Btu/ ft*, but the view 
factor has been decreased to 0, 70. On the basis of those values, much wisdom 
and knowledge wo»ild bo required to pro|>crly assess the effective radiator tem- 
perature. Note for high and a view factor of 0.70, the effective tempera- 

ture can bo outside the temperature range of the Inlet and outlet temperature. 
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APPENDIX A 

RADIATOR SIZING PROGRAM 


This computer program Is an equation solving procedure for one of two 
equations* The first is statement 49* which applies whenever the absorbed flux 
cannot be assumed zero* The second is statement 60* which applies when the 
absorbed flux is zero* Both equations are reported in Reference 2* 

Statements 3 through 27 determine the mass flow rate and resulting flow 
characteristics. 

The classical relationsldp between friction factor and Reynolds number 
is well known as the Moody diagram. It can readily be found in reference books 
on fluid flow* For laminar flow the Hagan- Poiseuille equation is transformed 
into the more manageable form shown on lino 12 whereas those values needed 
in the transition and turbulent regions (RE > 2100) require an iterative process* 
This is readily apparent from the Colebrook equation: 


= -0. 86 In 


^c/D 2* 51 \ 

^3.7 ^ RnTT/ * 


For a first approximation the right hand term containing the friction factor is 
Ignored and the resulting friction fat*tor is used for the next approximation. 
From this point a convergent routine is employed. These are illustrated in the 
program in lines 14 through 20* The pipe diameter is an input value and the 
roughness height is built into the program. Reynolds number is calculated in 
the classical manner from input values. The Darcy-Wcisbach equation is used 
to find the pressure difference which is used in the power equation. 

The following are definitions of the inpHit statement 2: 

TFl — Fluid Inlet Temperature, *F 
TF2 — Fluid Exit Tem|>eraturc, *F 
TWl — Wall Temperature at Inlet* *F 
TW2 — Wall Tcm|)craturc at Exit, *F 




f 

I ; 1 T 1 


.i - 4 



• ^ 



Q — Thermal Load» Btu/ sec 
CP — Heat Capacity of Fluid* Btu/ lb-* F 
imO — Mass Density of Fluid* Ibm/ft 
XMU — Viscosity of Fluid* Ibm/ft-hr 

XK - Conductivity of Fluid* Btu/hr-ft-*F 

X — Fluid Thickness (Duct Thickness) * In* 

E . Surface Emlsslvlty 

LD - Radiator Length Perpendicular to Flow, ft 

VP — View Factor 

CA — Absorbed Flux* Btu/hr-ft*-*F 

C — Case Identification Number 


The wall temperatures. TWl and TW2* are 
previous daU or other calculations. The program 
XW, statement 53 based on area. 


evaluated on the basis of 
computes the radiator length* 


B..lcny, the program oqoallons accouat for the temperature gr^leot 
Basically* p __*ancular duct necessitates a zero gradient 

*” the^^ Element 41 calculate, the flu effectltmuea. 

used should be no less than this optimum value. 
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OKIT n.ry FfkCB NATTS/FT 


tFlCt?} 



PRINT 19tTF 





8 ^ 

FCR»A^ riCTtlFM 

’CTAI 'IfW rpfcC* MATTS 


•riCtJi 

119* 


PRINT IttCITO 




119* 

19 

FOR«A^ lirTtSFlf 

Ftt'ir ’fFP. STFF, F 


•PlCtJl 

tlSt 


PRTN^ ITfPR 






f 

t 




r 


i i 

!• 



V 

t 


i 


J"i:TL U_jJ4^j J-J -i-LiJ . ij_jjjj :i i : . j 


18 



original PAGB to 
or POOR quautv 


IT 

rori<A? nrir»j<;H 
rUTMT 

»*AAAOt AU“Pt'’ 


•ric»3i 

It 

rOP*M 

^pcciri' ArrA 

rT?/p« 


10 

rOPPAT tlCY»2'*> 
30»0P 

rctrvAtrkT fa'“. 

TCPP. r 

tric.31 

?r 

ropHA’ f:cy*3^M 
ppi'it ntXMooT 

ppc?siJ»*r ?pcp 

rsi 

•ric«3i 

^l 

ropPAT acy*3fM 
rpinT 

PASt riry patt 

IC/HP 


?? 

reitPAT iicy#3*'K 
r»i«iT 2 S»cr 

PAS? vricciTT 

is/NP ni 

tric.3i 


rOPPAT ltCX*3CI 

ruPfAC riN err 

,tPTP. men 

ffSC.31 


ICO rCIt^AT IIHII 
1 VrOH*AT I I 
CO TO 111 
000 5T0» 

CNO 


Ml • I III lie Mill III Ml M I I II II Ml II Ml Ml I MMMMMM 

1 1 1 1 1 II I II III I II I II 1 1 1 1 1 1 1 1 1 1 1 1 II II 1 1 1 1 null III inn nn in iiniiniiin 
litUlitnUltlltJIt li l7JV7t77 7 7 2 UtU7VU 7 7 777777777777 77'7K7777^^ 7 7 7 7 
jj I ji 31 j 15 n 11 n n n n n n i >ii ii J J? n niminin 

41 1 44 4 444 1444 4444 4 I 4 4 144 1 t4<44«4l 44 4 4 I< 1(1 < 4 ( ( 4 4 4II444444444444II 1 : 14 4 114444 1 4 

ssisssnss SiS 5 SSiSS 5 S>>n 5 iSJSiS‘ 5'5 5 5 i 5 S 5 S‘jSniSli 5 J$iS 5551 in 5 SnjM$lJiS$S 

I 7 7 7 I 7 7 7 i/)!i77ll7ni/l/l/l t ) 1 1 7 7 7 i 7 1 7 7 7 1 1 1 1 7 7 7 7 7 1 1 7 7 7 7 7 7 7 7 7 7 7 1 7 7 7 ! 1 1 7 7 7 7 7 7 7 7 7 

II I II II I $! I! II I M if ti n { t!l ill 4 II M Mllllllim 

»*iiii4i miiii'''*4i» jj 14 ni*”'* 4»44 *'>iii.4>sss5?'Hiii^4^imi** '•Ji 

g I ‘ 4 /, », ; , 4» ■ ». • . . ; ; • 4'- t »m > i.i i »■.?». %r . • • •! 


CX3J3J M il I -I M JJ l i'l i 




I 


1 


I 


I 


INPUTS 


CASC I*0* N0«****»**«*** II»C1C 


rtUTO INLET TCNP, F 97.C0C 
FLUID EXIT TEHP. F «0.000 
RAO. NALL TENP. 9 INLET F 95.PCC 
RAD. MALL TENP. D OUTLET F St.OOO 
RAD. THERMAL LOAD BTU/SEC lA.IIO 
fluid heat CAPCITT 8TU/LB F .210 
FLOTO OENSITT L8/FT3 87.CCC 
FLUID VTSOCITY LB/FT HR .90J 
'^LOTO CONDUCTIVITY BTL'/HR FT F ,C9C 
FLUID EOU. THICKNESS IN .100 
RAO. SURFACE ENISSIVITY ,9CC 
RAO. LENGTH PER. TO FLOM FT 12.000 
VIFM factor l.CCO 
total ENVIR. absorbed BTU/HR FT2 15.000 


OUTPUT 


REOD. REJECTION AREA FT 2 

R7B.2T 

RADIATOR SIZE FT 2 

739. IR 

RADIATOR LENGTH FT 

19.93 

REYNOLDS NUMBER 

751. 038 

FILM COEFF. BTU/HR FT2 F 

35 .088 

UNTT FLOW POWER WATTS/FT 

.C02 

TOTAL FLOM POWER MATTS 

.001 

FLUID TEHP. OIFF. F 

57. roc 

PRANOL NUMBER 

7.200 

SPECIFIC AREA FT2/KW 

31.905 

EQUIVALENT RAO. TEHP. F 

53.237 

PRESSURE DROP PSI 

.CI6 

NASS FLOW RATE LB/HR 

0082.297 

NASS VELOCITY LB/NR FT2 

0C822.988 

SURFACE FIN EFF.tHIN. WEIGHT 

•502 
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APPENDIX B 

RADIAKIR TEMPERATURE DISTRIBUTION 

In a differential radiator length, dL, the radiator temperature will change 
dT in accordance with the following energy balance: 

jc^T^ - cj Lj dL = -ACp dT . (B-D 

This equation is presented in Reference 2 by Mackey. In this form it assumes 
the fluid temperature is the same as the radiator wall temperature. If the 
environmei al factor, C , can be assumed zero, the equation can be readily 

integitited: 



in thia form, at a point on the iradiator having temperature, T, the iradiator 
length must be L. 

Equation (B-2) can be combined with equation (12) to )rlold 



In this form, L/L is the decimal value of the total length, L . The radiator 
o 

Umperature, T, corresponds to the decimal length, 
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